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ABSTRACT 


We present a practical importance-sampling scheme for the Student- 
T distribution of visible normals by representing the Student-T NDF 
as a superposition of Beckmann NDFs. Additionally, we derive a 
new form of delta tracking to evaluate and sample exact BSDFs with 
general full-sphere NDFs. These tools permit efficient computation 
of benchmark BSDF values for the multiple scattering from general 
(including porous) rough surfaces. 
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1 INTRODUCTION 


The normal distribution function (NDF) of microfacet theory strongly 
influences the appearance of rough surfaces. Beckmann and GGX 
NDFs are most commonly used in rendering, but do not describe 
all surfaces, and more general distributions are required in some 
cases. The two-parameter Student-T NDF was introduced to include 
Beckmann and GGX within a broader family [Ribardière et al. 2017], 
but has not been compared to measured data. Further, Student-T 
is limited to describing heightfields, as are other shape-invariant 
NDFs like the Bessel-K NDF [Bahar and Fitzwater 1983]. To model 
very rough and porous surfaces, an NDF that extends to the full 
sphere is required [Dupuy et al. 2016]. 

Despite the motivation for more general parametric NDFs, in 
practice only Beckmann and GGX enjoy exact visible-distribution- 
of-normals (vNDF) sampling procedures that permit low-variance 
rendering in a path tracer with accurate multiple scattering via 
random walks. To extend vNDF sampling and random walk multiple 
scattering to more general NDFs, we present a number of new BSDF 
tools. Specifically, in this talk we: 

(1) Show that Student-T very accurately matches the behaviour 

of frosted and etched glass (Figure 1) 
(2) Derive a vVNDF sampling scheme for Student-T 
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(a) Frosted Glass 


(b) Etched Glass 


Figure 1: The Student-T NDF produces very close matches 
to the etched and frosted samples measured previously by 
[Walter et al. 2007] for y = 3, outperforming both Beckmann 
and GGX fits. 


(3) Derive a new form of delta tracking that can produce bias- 
free reference solutions for surfaces with any bounded NDF 

(4) Release a new codebase with validated benchmark implemen- 
tations of the above methods, behind a new BSDF interface 
that is easily extendable to new NDFs and permits any BSDF 
on the microfacets (allowing bi-scale roughness). 


2 THE STUDENT-T NDF 


2.1 New Physical Motivation 


The Student-T NDF was primarily motivated by convenient mathe- 
matical properties that are required for efficient use (such as closed- 
form shadowing/masking factors). To the best of our knowledge, 
ST BSDFs have not been compared to measured data in order to 
determine if ST is indeed a good physically-motivated choice for 
interpolating between Beckmann and GGX. We revisit the original 
data of transmission through rough glass by [Walter et al. 2007] and 
show that the ST NDF is indeed a very strong choice, matching the 
etched and frosted samples very closely. We highlight the accuracy 
of these two fits in Figure 1 and refer the reader to the supplemental 
text for additional details. 


2.2 Sampling the Student-T Distribution of 
Visible Normals 


With notation u = cos 0m, the Beckmann NDF of roughness a is 


ce 


e ‘ae 
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where ©(u) is Heaviside’s function. The Student-T (ST) NDF with 
shape parameter y > 3/2 is 


(y -= 1)? (æy -1)+4- 1)” 
. ou). (2) 


DST (u, ay)= a 
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Figure 2: Illustration of our stochastic null-scattering formu- 
lation of non-heightfield Smith scattering. 


The key insight to our method is to note that the ST NDF can be 
expressed as a gamma-superposition of Beckmann NDFs: 


co y- —1 
D"(u, ay) = f fae u, a] i) dm. (3) 
0 m 


r(y-1) 
Here m > 0 indexes the superposition, and for each m the rough- 
ness of the corresponding Beckmann distribution in the mixture is 
ayy- 1I)/m. 
This superposition can only be used directly for normal incidence 
(u = 1). In the supplemental material, we show how a related 


=m 2 


DE 


distribution f (m, uj) can be used to sample superposition roughness 
m given incident cosine uj such that sampling Beckmann vNDFs 
with m results in exact VNDF sampling of ST. New renders of rough 
gold and diffuse surfaces with multiple-scattering ST BRDFs are 
included in the supplemental material. 


3 ANULL-SCATTERING FORMULATION OF 
ROUGH SURFACE SCATTERING 


We also derive a new Monte Carlo procedure for stochastically 
evaluating and sampling Smith microfacet BRDFs. Our method lifts 
the heightfield requirement of classical theory (that all microfacet 
normals lie in the upward hemisphere), allowing porous and purely 
volumetric materials to be described together with heightfields 
in a single unified formalism. Our approach takes as input any 
bounded NDF on the full sphere. Random walks in the volume are 
then sampled after introducing ficticious/null microflakes over the 
space of microfacet orientations such that the new NDF is uniform 
over the sphere, resulting in uniform cross sections and uniform 
distributions of visible normals. The complex per-NDF derivation 
of distance-sampling and visible-normal sampling routines of prior 
work is thereby avoided. The tradeoff is additional cost that is re- 
quired to perform a rejection procedure analogous to null-collisions 
for distance sampling in heterogeneous volumes. We show that 
this approach exactly includes prior height-field Smith models as 
special cases, but note that the efficiency greatly suffers for low- 
roughness materials. We anticipate that the approach will be useful 
for LOD compression of geometric complexity and for modeling 
very rough / porous / granular / particulate BRDFs (to benchmark 
faster approximations or train neural models). 
In summary our method works by: 
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e Given a bounded NDF D(wm) 2 0, we select a constant 
majorant NDF (wm) > D(wm) over S?. 

e We begin a random walk at the boundary of a half space, 
along w = —oi. 

e We sample collisions in the medium using the majorant NDF 
D(@m), which reduces to a constant, allowing o; = 1. 

e While no actual particles are traced against, the following 
geometrical procedure illustrates the core sampling idea. 
For each sampled collision (real or ficticious), we consider 
a small spherical particle with the pair of NDF values (real 
and ficticious) mapped to its surface (Figure 2). We then 
consider a uniform field of rays along direction w intersecting 
this sphere, and sample from this set of rays uniformly at 
random (2D disk sampling). The intersection of this ray with 
the sphere produces a microfacet normal wm: the normal 
at the intersected location. This sampled collision is real 
with probability D(@m)/D(@m) and otherwise we ignore 
the collision and sample a new distance. 

e For each real collision, the microfacet BSDF is sampled using 
w and wm to determine the next direction. 

e Upon exiting the half space, the particle weight (which may 
deviate from 1 for absorbing microfacets) and direction w 
are returned. 

Stochastic evaluation of the BSDF performs the same random walk, 
but with NEE along wo at each real collision. We numerically es- 
timate the Smith Lambda function for any NDF using a Green’s 
function (the Lambda function for the Dirac delta NDF). See the 
supplemental pdf and code for additional details. 


4 CODE AND RESULTS 


In the supplemental material we provide a validated reference C++ 
implementation of the above BSDF tools together with a number 
of images showing novel BSDFs, such as ST rough metals with 
multiple scattering, biscale roughness and porous metals with a 
vMF NDF. The code is implemented behind a new compact interface 
that allows passing any BSDF and NDF into a common Microsurface 
class, resulting in a new BSDF that can be recursively applied to 
another microsurface for biscale roughness. 
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